Absorption losses in periodic arrays of thin metallic wires 
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We analyze the transmission and reflection of the electromagnetic wave calculated from transfer 
matrix simulations of periodic arrangements of thin metallic wires. The effective permittivity and 
the absorption is determined. Their dependence on the wire thickness and the conductance of the 
metallic wires is studied. The cutoff frequency or effective plasma frequency is obtained and is 
compared with analytical predictions. It is shown that the periodic arrangement of wires exhibits 
a frequency region in which the real part of the permittivity is negative while its imaginary part is 
very small. This behavior is seen for wires with thickness as small as 17 /im with a lattice constant 
of 3.33 mm. 



Rapidly increasing interest in left-handed metamateri- 
als (LHM) (for a recent review, seei) raised some interest- 
ing questions about the electromagnetic (EM) properties 
of composites which contain very thin metaUic compo- 
nents. The most simple example of such a composite is 
a periodic array of thin metallic wires. Pendry el ali^ 
predicted that such system behaves as a high pass filter 
with an effective permittivity 
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In Eq. n fp is the effective plasma frequency, or cutoff 
frequency^, and 7 is the damping factor. Various theo- 
retical formulas were derived for the dependence of the 
plasma frequency on the lattice period a and the wire 
radius r. Pendry el al^ obtained that 
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Shalaev and Sarychev obtained that 
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In Eqs. I2I4I cught is the velocity of light in vacuum. 

Periodic arrangements of thin metallic wires are 
used as a negative-e mediumi^ in the left-handed 
structurea^iifliii. It is therefore important to understand 
how the electromagnetic response - not only the effec- 
tive plasma frequency, but also the factor 7 - depends 
on the structural parameters of the wire system. Re- 
cently, Ponizhovskaya el alP claimed that for small wire 
radius, the absorption in the wire system is so large that 
the transmission losses do not allow any propagation 
of EM wave in left-handed structure. Very low trans- 
mission, measured in the original experiments on the 



LHMpiiS seemed to agree with their pessimistic conclu- 
sion. However, it is not clear why the transmission was 
so low in the original experiments. Recent experimental 
measurements-'^^'''^^- established that the transmission of 
LHM could be as good as in the right-handed systems. 

Our aim in this Letter is to study numerically how 
the effective permittivity of the periodic arrangement of 
metallic wires depends on the wire radius and on the 
conductance of the wires. We present results for the real 
(e^ff) and imaginary (e"ff) part of the effective permittiv- 
ity of the wire medium, estimate the transmission losses 
and the plasma frequency and compare our results with 
the analytical formulas given in Eqns. I1I4I 

In our numerical simulations we use the transfer ma- 
trix method (TMM). Details of the method are given 
elsewhere-^^. Here we only point out the main advantage 
of the TMM, namely that it gives directly the transmis- 
sion t, reflection r and absorption A = 1 ~ |ip — |rp of 
the EM plane wave passing through the system. Con- 
trary to this, in the FDTD method, which was used in 
Ref^, one obtains t and r form the time development of 
the wave packet, which is much more complicated and 
probably also less accurate. To be able to obtain eos, 
one also needs the phase of r and t, in addition to their 
amplitudes. This is also easily achieved in the TMM. 

From the obtained data for the transmission and re- 
flection, we calculate the effective permittivity of the 
system.^* The refraction index is given by the equation 
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Since we do not expect any magnetic response, we fixed 
the value of the permeability to be /i = 1. The permit- 
tivity is then found as eoff = n^. 

The way we discrete the space in the TMM might con- 
trol the accuracy of our results. To test how discretiza- 
tion influences our results, we repeated the numerical 
simulation for different discretizations. The wire is repre- 
sented as a rectangular with square cross section 2r x 2r, 
r being the "wire radius". The obtained results for the 
effective permittivity CeS (both real and imaginary part) 
are almost independent of the discretization procedure. 



2 



Figure 1 shows how the effective permittivity depends 
on the wire radius. We analyzed four different wire ar- 
rays with period a = 3.33 mm. For all of them the real 
part of the effective permittivity is negative and could be 
fitted by Eq. (QJ. This enables us to obtain easily the 
plasma frequency. Only for the smallest wire thickness 
studied (17 x 17^m) one gets a relatively large e"ff for 
small frequencies. In this region relation (^) is not valid. 
Nevertheless, for frequencies larger than 5 GHz, e"jj is 
small and e'^g negative. 

The right lower panel of fig. 1 shows data for wires with 
the cross section of 17 x 300/j,m. These parameters were 
used in the experiment of Shelby et ali^. We again see 
that the formula given by Eq. qualitatively agree with 
our data. Thus, there is no doubt that this array of wires 
really produces a medium with negative e^g, which then 
can be used in the creation of the left handed systems. 

Figure 2 compares our data for plasma frequency with 
the analytical formulas given by Eqns. I2I4I Accepting 
some uncertainty in the estimation of the plasma fre- 
quency from the numerical data, we can conclude that 
for thin wires our data agree with the theoretical formula 
(O of Shalaev and Sarychev4. For thicker wires, our re- 
sults are in agreement with the formula Q of Maslovski 
et al^ 

We also study the role of the conductance of the metal- 
lic wires. In the simulations shown in fig. 1, we consider 
the metallic permittivity to be — (—3 -f 588 i) x 10^. 
We are aware that this value of em is smaller than the 
permittivity of realistic metallic wires: for instance for 
copper Cm ~ 5 x 10^ z, as follows from the relation be- 
tween the permittivity and the conductancei^ (the con- 
ductivity of copper is a « 5.9 x l{f {U,m)~^) . Our data in 
fig. 1 therefore underestimate losses, because transmis- 
sion losses are smaller for higher values of the metallic 
permittivitjii^. This is clearly shown in fig. 3 where we 
present the ratio k = \e'^g/e'^g\ vs frequency for two sys- 
tems which differ only in the value of the imaginary part 
of the metallic permittivity. Fig. 3 shows also the fre- 
quency dependence of the absorption as obtained from 
the numerical simulations. The absorption also exhibits 
a maximum in the neighbor of the plasma frequency. 

Fig. 1 also confirms that e'^^ increases when the wire 



radius decreases. For instance, the parameter 7 is only 
0.003 GHz for r = 100/im. As we will see below, (fig. 4) 7 
increases up to 1.2 GHz when the wire radius decreases to 
15/im. Nevertheless, even for very thin wires, losses are 
much less than what was claimed in Rcf . ^ . As it is shown 
in fig 1, also an array of wires with thickness 17 x 17/im 
creates a negative-e medium. As this is in strong contrast 
with the results of Ref.^, we decided to study exactly the 
same system as that o&. Results of our simulations are 
shown in fig. 4. Although such systems are not used in 
experimental arrangements of the LHM, our results give a 
comparison between two different numerical treatments. 
Our data again clearly show that e^g is negative for / < 
fp. As it is shown in the inset, the transmission losses 
are also small. 

We believe that the present data are more accurate, 
than those published previously in^, not only because 
they agree with the theoretical analysis but also because 
the TMM gives the reflection and its phase straightfor- 
ward. This is important because the main difference be- 
tween our results and those of& seems to be in the estima- 
tion of reflection R = |rp. When comparing our data for 
absorption, given in the inset of fig. 4 with those given 
in fig 2b of Refi^, we see that our absorption is much less 
than that estimated in Refi^. 

In conclusion, we analyzed numerically the transmis- 
sion properties of a periodic arrangements of thin metallic 
wires. From the transmission and reflection data we cal- 
culate the effective permittivity and plasma frequency, 
which agrees qualitatively with theoretical predictions. 
Both the effective permittivity and the absorption data 
confirm that the array of thin metallic wires, used in the 
recent experiments on the left-handed metamaterials, in- 
deed behaves as a negative permittivity medium with low 
losses. 
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FIG. 1: Effective permittivity as a function of frequency for 
various shapes of the metallic wires. The lattice period in 
all cases is a = 3.33 mm. We used the metallic permittivity 
Em = (-3 + 588 i) X 10^ 
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FIG. 2: Plasma frequency as a function of the wire radius. 
The lattice constant is a = 5 mm. Solid, dashed and dot- 
dashed line is the result of Pendry et al. @, Shalaev and 
Sarychev Q and Maslovski et al. respectively. Dot line 
is a fit of our data to the function /p = ao / ^J\n{a\/r) with 
parameters ao = 20.9 and a\ — 0.84. 
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FIG. 3: Ratio k — |e"ff /ecff | for a lattice of wires with radius 
50/im. Metallic permittivity is em = (—3 + 588 i) x 10'^ (open 
circles) and = (-3 + 5 880 i) x 10^ (full circles). Dashed 
(solid) line is absorption for the corresponding systems ob- 
tained numerically by the TMM. These numerical results con- 
firm that losses are smaller for higher metallic permittivity 
and that the value of the plasma frequency, estimated ap- 
proximately from the position of the maximum of k, does not 
depend on the value of the metallic permittivity. 




FIG. 4: Effective permittivity for a lattice of thin metallic 
wires. The wire radius is r = 15pim and the lattice constant 
is o = 5 mm. The metallic permittivity is em = —2000 + 
10^ i. Two different discretizations are used with mesh sizes of 
30/im (open symbols) and 15pim (full symbols) . The solid and 
dashed lines are fit to Eq. 0with fp — 11.1 GHz and 7 = 1.2 
GHz. The length of the system was up to 60 unit lengths 
for open symbols and 10 unit lengths for full symbols. Inset 
shows the numerically calculated absorption as a function of 
frequency. 



